x=1:1:10;
y=x.^2;
plot(x,y,'+')